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Abstract 

The response of mesophyll conductance to C0 2 (g m ) to environmental variation is a challenging parameter to 
measure with current methods. The 'variable J' technique, used in the majority of studies of g m , assumes a one-to- 
one relationship between photosystem II (PSII) fluorescence and photosynthesis under non-photorespiratory 
conditions. When calibrating this relationship for Populus trichocarpa, it was found that calibration relationships 
produced using variation in light and C0 2 were not equivalent, and in all cases the relationships were non- 
linear—something not accounted for in previous studies. Detailed analyses were performed of whether different 
calibration procedures affect the observed g m response to C0 2 . Past linear and assumed calibration methods 
resulted in systematic biases in the fluorescence estimates of electron transport. A sensitivity analysis on modelled 
data (where g m was held constant) demonstrated that biases in the estimation of electron transport as small as 2% 
(—0.5 |mmol m~ 2 s~ 1 ) resulted in apparent changes in the relationship of g m to C0 2 of similar shape and magnitude to 
those observed with past calibration techniques. This sensitivity to biases introduced during calibrations leads to 
results where g m artefactually decreases with C0 2 , assuming that g m is constant; if g m responds to C0 2 , then biases 
associated with past calibration methods would lead to overestimates of the slope of the relationship. Non-linear 
calibrations were evaluated; these removed the bias present in past calibrations, but the method remained sensitive 
to measurement errors. Thus measurement errors, calibration non-linearities leading to bias, and the sensitivity of 
variable J g m hinders its use under conditions of varying C0 2 or light. 

Key words: Chlorophyll fluorescence, curve fitting, electron transport rate, g m , mesophyll conductance to C0 2 , Populus 
trichocarpa, variable J technique. 



Introduction 

Mesophyll conductance (g m ) is the conductance of C0 2 Whether g m is a constitutive or dynamic characteristic of 

from the intercellular airspaces to Rubisco, a largely liquid a leaf is fundamental to our understanding of plant 

pathway through the cell wall and three membranes, responses to the environment. As g m may represent up to 



Abbreviations: a, leaf absorptance; (3, fraction of quanta absorbed by PSII; 0 R3 n, quantum efficiency of PSII; <Dco 2 , quantum efficiency of gas exchange; r*, photo- 
compensation point; A, net photosynthetic rate; A c , Rubisco-limited photosynthetic rate; A Jt RuBP regeneration-limited photosynthetic rate; /4 TPU , triose phosphate 
utilization-limited photosynthetic rate; C c , chloroplastic C0 2 concentration; C„ intercellular C0 2 concentration; C*, apparent photo-compensation point; EDO, 
exhaustive dual optimization procedure; g m , mesophyll conductance to C0 2 ; K c , Rubisco Michaelis-Menten constant for carboxylation; K Q , Rubisco Michaelis-Menten 
constant for photorespiration; J, electroxn transport rate; J A +R d » rate of J needed to account for measured A+R d ; J cal> , calibrated fluorescence-derived J; J raw , 
uncalibrated fluorescence-derived J; J to tai> modelled or measured total J to carboxylation and photorespiration; O, oxygen mole fraction; PPFD, photosynthetic photon 
flux density; R d , mitochondrial respiration in the light; S c/0 , relative specificity of Rubisco; T h leaf temperature; V c , rate of carboxylation; V CiC , rate of carboxylation limited 
by Rubisco; V c j, rate of carboxylation limited by RuBP regeneration; \Z Cim j n , minimum of V c>c and V cj ; V cmax , maximum rate of Rubisco carboxylation; V Q , rate of 
photorespiration; V alh alternative J in C0 2 equivalents; VPD, vapour pressure deficit. 
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40% of the C0 2 diffusional limitation on photosynthesis 
(Warren, 2008), dynamic variation in g m would offer a major 
avenue for photosynthetic regulation, comparable with that 
of the stomata. The most commonly used technique to 
measure g m , the variable J method, consistently demon- 
strates a large reduction in g m with increasing C0 2 (Flexas 
et al, 2007). However, the size and presence of the response 
of g m to C0 2 varies between studies using a variety of 
methods (Flexas et al , 2007; Tazoe et al , 2009; Vrabl et al , 
2009). For example, a less steep response of g m to C0 2 was 
found for Nicotiana tabacum using the independent carbon 
isotope method relative to the variable J method (Flexas 
et al, 2007). In a separate experiment, Arabidopsis thaliana 
and N. tabacum were reported to reduce g m by —85% and 
65%, respectively, when C0 2 changed from 200 umol mol -1 
to 1000 umol mol -1 at 21% 0 2 and measured using the 
variable J method (Flexas et al, 2007). In a second 
investigation, the carbon isotope method resulted in only 
a 10% reduction and a 5% increase for the same species, 
respectively, when measured across the same range of C0 2 
mole factions; measurements at 2% 0 2 showed reductions 
of 26% and 40% (Tazoe et al, 2011). The widely used curve- 
fitting techniques for estimating g m explicitly assume 
a constant g m across the range of C0 2 used to generate 
C0 2 response curves (Ethier et al, 2006; Warren, 2006; 
Sharkey et al, 2007; Gu et al, 2010). To date the 
underlying mechanisms determining g m and the reason for 
the different results between methods remain unresolved. 

The 'variable f technique encompasses a group of 
methods that estimate g m , chloroplastic C0 2 concentration 
(C c ), and the rate of oxygenation or photorespiration (V 0 ) 
from combined fluorescence and gas exchange data. Meso- 
phyll conductance to C0 2 is calculated as the ratio of net 
photosynthetic C0 2 flux (A) to the difference in C0 2 
concentration between the intercellular airspaces (Q) and 
the chloroplast (C c ). C c is related to the ratio of carboxyl- 
ation to photorespiration at Rubisco, and photorespiration 
is then proportional to the difference between fluorescence- 
derived estimates of the total electron transport rate and the 
rate of electron use by carboxylation estimated from gas 
exchange. This derivation is described in detail in the 
Materials and methods and reviewed by Warren (2006) and 
by Pons et al (2009). 

Fluorescence estimates of total electron transport are 
derived from the work of Genty et al (1989), who 
established that under non-photorespiratory conditions 
(low oxygen and high C0 2 ), a linear relationship exists 
between the quantum yield of fluorescence (® PSII ) and 
measured quantum efficiency of rates of C0 2 fixation 
(OC0 2 ). This proportionality has subsequently been used 
to provide an estimate of electron transport rates: in the 
absence of alternative electron sinks, the relationship 
between carboxylation estimates of linear electron flow 
and fluorescence estimates of electron transport should be 
one-to-one under non-photorespiratory conditions. In 
practice, this relationship deviates from one-to-one due to 
interspecific variation in the values of standard constants 
such as leaf absorptance (Baker, 2008), and measurement 



of the relationship under non-saturating C0 2 where 
significant alternative electron transport sinks may be 
present. However, for simplicity, it is often assumed that 
standard constants are accurate and do not vary during 
experiments. 

An alternative approach is to conduct pre-experimental 
calibrations to provide estimates of electron transport from 
photosystem II (PSII) fluorescence (Lawlor and Tezara, 
2009). While empirical calibration has the potential to 
improve estimates of electron transport, it can also in- 
troduce systematic errors (biases) in the calculation of the 
total electron transport rate, and thus g m . The impact of 
calibration issues, such as non-linearity, on the calculation 
of g m has not been thoroughly assessed. 

The present study examined whether the differences 
between two common methods used for measuring g m is 
the result of biases in the calibration of the variable 
J method. However, the challenges inherent in the variable 
J method have long been recognized (Harley et al, 1992), 
with the Harley criterion providing an indication of how 
sensitive the g m values are to errors when using this method 
(Harley et al, 1992). The original sensitivity analyses 
of Harley et al (1992) demonstrated that the relationship 
of g m to C0 2 was sensitive to errors in the values of 
mitochondrial respiration, the photo-compensation point, 
and the fluorescence estimate of the total electron transport 
rate. However, this analysis was not extended to a broad 
range of Qs, as subsequent studies do, and the sensitivity of 
the g m response to C0 2 to errors has not been compared 
with the size of biases present in the calibration procedure. 

Therefore, the goal of this study was to understand the 
conditions under which g m can be accurately measured 
using the variable J technique. Fluorescence with gas 
exchange measurements is calibrated using classical meth- 
ods for the widely used genome model plant poplar 
(Populus trichocarpa Torr. & Gray). Consistent with the 
original literature, significant variation in the calibration 
relationship was found, such that there is the potential for 
systematic error when calibrations are applied to a broad 
range of environmental conditions. Photosynthetic models 
are then used to assess the effects of biases on the response 
of g m to C0 2 . Finally, new calibration techniques by which 
these biases may be reduced when estimating a single value 
of g m for a leaf, or comparing species, are suggested. 
However, it is demonstrated that the variable J method 
should be used with caution when measuring the response 
of g m to C0 2 and light, as any bias in the estimation of 
electron transport rates results in changes in the relationship 
of g m with C0 2 . 



Materials and methods 

Plant material and growing conditions 

Poplar plants were propagated from cuttings and grown in 
environmentally controlled growth chambers. Metal-halide and 
high pressure sodium lighting (400 umol m -2 s _1 ) was provided 
for 14 h per day. Temperatures in the chambers were maintained 
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between 20 °C and 24 °C, and humidity was kept at 70%. The 
cuttings were placed in 3785 cm 3 pots in Farfard 3B potting soil 
which included Osmocote Plus slow release fertilizer as per the 
manufacturer's instructions (15/9/12/1 N/P/K/Mg plus trace ele- 
ments: S, B, Cu, Fe, Mn, Mo, and Zn; Scotts Company, OH, 
USA). The pots were watered daily and fertilized weekly with 
Peters Excel All Purpose soluble fertilizer (21/5/20 N/P/K plus 
trace elements: B, Cu, Fe, Mn, Mo, and Zn; Scotts Company). 
Plants were measured after 4-9 months of growth in March-May 
of 2010 (all experiments), and a second set of plants in November 
2010 (extra C0 2 response curves). 

Gas exchange and fluorescence measurements 

Gas exchange and fluorescence measurements were done on young 
fully expanded leaves using a 2 cm 2 LI-COR LI-6400 fluorescence 
chamber and gas exchange system (LI-COR, Lincoln, NE, USA). 
Plants were allowed to acclimate to the gas exchange system in 
a laboratory growth chamber for >30 min, until stomatal conduc- 
tance was stable. Unless otherwise noted, general measurement 
conditions were as follows: photosynthetic photon flux density 
(PPFD), 400 umol m -2 s _1 with no blue light component (Loreto 
et aL, 2009); reference C0 2 , 400 umol mol -1 ; T h 24.9±0.8 °C; 
vapour pressure deficit (VPD), 1.47 ±0.49 kPa; flow, 150 umol s" 1 . 
Single flash fluorescence measurement settings were used and 
adjusted according to the optimal values obtained from the flash and 
measuring intensity procedures in the LI-6400 manual (Anon, 2004). 
All measurements were corrected for leaks using empirically de- 
termined leak corrections for dry poplar leaves under measurement 
conditions (C02S correct ed=C02S uncorrect e d +9.868x 10~ 7 xCO2R 2 + 
5.291xl0~ 4 xCO2R-1.008). 

The Laisk method was used to measure non-photorespiratory 
respiration in the light (R^) and the apparent photo-compensation 
point (Q*) (Warren, 2006). In this method, the y and x value of 
the average intersection of three C0 2 responses are taken as 
and Q*. The C0 2 responses were measured at reference C0 2 
concentrations of 150, 100, 75, and 50 umol mol -1 , at three light 
levels 400, 175, and 75 umol m -2 s -1 . Six replicate sets of 
measurements were made providing mean (±SD) values of 
0.42±0.21 umol m -2 s -1 for R d , and 36.0±1.9 umol mol" 1 for 
Q*. In theory, the Cf values should be increased by RJg m to 
obtain an estimate of the photocompensation point T* (von 
Caemmerer, 2000), but as no independent value for g m was 
available, the transformation was not performed and T* was taken 
to be equal to Q*. A sensitivity analysis, described below, 
confirmed that minor variation in the value for T* did not greatly 
affect the values of g m (Table 1). 

Low oxygen, non-photorespiratory conditions, were obtained by 
mixing air with nitrogen gas using a Wosthoff gas mixer to achieve 
a 1% 0 2 content. This was tested using C0 2 drawdown as an 
indicator of the 5% mixing ratio necessary to produce 1% 0 2 from 
21% 0 2 air, and further verified using an Ocean Optics USB4000- 
FL-450 Fiber Optic Spectrophotometer oxygen probe. Initial 
experiments demonstrated that 1% 0 2 was the highest oxygen 
concentration at which the slope of the variable J calibration 
relationship did not change with successive dilution of 0 2 at 
ambient C0 2 , indicating that further reductions in 0 2 would not 
further inhibit photorespiration. 

Calibration relationships and C0 2 responses 

Light or C0 2 response curves were measured under non-photo- 
respiratory conditions to calibrate the relationship between 
fluorescence-derived electron transport rates (/ raw ) an d photosyn- 
thesis (J a +R d ). Nine light response curves were measured at 
ambient C0 2 (400 umol mol -1 ) and 1% 0 2 , starting at a PPFD of 
2000 umol m -2 s -1 and reducing PPFD to 1500, 1000, 800, 600, 
500, 400, 350, 300, 250, 200, 150, and 100 umol m -2 s -1 with 
4-6 min intervals between measurements. Five C0 2 response 
curves were measured at 400 umol m -2 s -1 PPFD and three at 



Table 1 . Sensitivity analysis for the calibration and calculation of 
g m at ambient C0 2 (400 iimol moi~ 1 ) and 400 }xmol m~ 2 s~ 1 
PPFD 

All percentage increases were expressed relative to the mean for 
scenario 1a which used standard calibration parameters (oc=0.85, 
P=0.5). Scenarios 1b-1e represent g m values calculated using the 
standard calibration parameters and measured variation in R d and 
r*. Scenario 2 represents the shift in g m due to measuring, and not 
assuming oc=0.85. Scenarios 3a and 3b represent g m values 
calculated using calibrations estimated from different plots and 
regression fits to light response curves. Scenario 4 represents g m 
values calculated using variation in C0 2 to generate a linear- 
sigmoidal calibration curve. Values represent means and standard 
errors for six replicates. 



Calibration Parameter Apparent g m % 

type varied (mol m~ 2 s~ 1 ) increase 

with 

calibration 



No calibration standard parameters 



1a oc=0.85, p=0.5 


NA 


0.166±0.018 


0 


1boc=0.85, (3=0.5 


R d +95% CI 4 


3 0.173±0.020 


4.1 


1c a=0.85, (3=0.5 


/V95% CI 


0.160±0.017 


-3.7 


1d 01=0.85, (3=0.5 


T*+95% CI 


0.181 ±0.021 


9.0 


1e a=0.85, (3=0.5 


r*-95% CI 


0.154±0.016 


-7.5 


Measured leaf absorptance 








2 oc=0.831 (measured), 


NA 


0.182±0.022 


9.4 


(3=0.5 








Light response at 1 % 0 2 and ambient C0 2 


calibration 




3a Linear fit 5 


NA 


0.341 ±0.087 


104.5 


3b Linear-sigmoidal 


NA 


0.306±0.079 


83.6 



C0 2 response at 1 % 0 2 and 400 \\n\o\ m 2 s 1 light calibration 
4 Linear-sigmoidal NA 0.199±0.017 19.8 

a R 6 and r* were estimated from six replicate sets of Laisk curves: 
in 1 b-1e the mean values were used plus or minus the 95% 
confidence interval of the mean for the replicates (R d , 0.42 ±0.1 7 umol 
m" 2 s" 1 and r*, 36.0±1 .5 umol mol" 1 ); b Linear fit to high light data 
from the efficiency plot (Fig. 1C) for data where O c02 is <0.05. 

1000 umol m -2 s -1 PPFD at 1% 0 2 . C0 2 was reduced from 
400 umol mol -1 to 100 umol mol -1 in decrements of 75 umol 
mol -1 , and after an 8 min re-acclimation at 400 umol mol -1 
increasing C0 2 to 600, 800, 1000, 1500, and 2050 umol mol -1 . 
Gas exchange was measured at each C0 2 concentration after the 
cuvette C0 2 concentration was stable for >120 s. A second similar 
series of C0 2 responses was measured after the first under 21% 0 2 . 
The leaf absorptance of 10 leaves was measured using a Taylor 
integrating sphere (LI-COR 1800-12). 



Estimation of 'variable J' g m 

Values for g m were estimated from the following standard 
formulae used in the 'variable f technique (Harley et aL, 1992; 
Valentini et aL, 1995; von Caemmerer, 2000), and using 
variants of the calibrations detailed below. Mesophyll conduc- 
tance to C0 2 is estimated as the ratio of the net photosynthetic 
rate (A) and the difference in C0 2 mole fraction from the 
intercellular airspaces (Q) and the chloroplastic sites of 
photosynthesis (C c ): 

As A and Q are provided by standard gas exchange measure- 
ments, estimation of C c remains as the difficult-to-measure un- 
known in this equation. C c is estimated assuming that Rubisco 
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specificity to 0 2 and C0 2 (S c / 0 ) remains constant, and that the 
ratio of the carboxylation rate (V c ) to the oxygenation rate (V 0 ) 
varies in direct proportion to the concentration of C0 2 at the site 
of carboxylation (C c ) in the chloroplast or the concentration of 
oxygen (O) which is assumed not to vary. Thus: 



where T* is the photo-compensation point (=0.5xOAS c / o ), mea- 
sured using the Laisk method. V c can be estimated as the sum of 
the measured A, a value for R d assumed to be constant and equal 
to that measured using the Laisk method, and half of V Q \ 

V c = A + R d + 0.5V o (3) 

V 0 is included in V c , as for every two oxygenations one C0 2 is 
released, leading to gross photosynthesis being underestimated by A. 
The total electron transport (/total) is tne sum °f tne reductant 
required for V c , V 0 , and any alternative electron transport sinks 
(^ait)- Under many conditions four electrons are used per C0 2 
molecule fixed (Baker, 2008), and it is known that two photo- 
respiratory cycles release one C0 2 , thus the rate of photorespiration 
can be estimated by rearranging this equation: 

V 0 = /total/4 - ( V c + K alt ) = ^(/ t otal/4 -{A + RJ+ V&)) (4) 

Here / total includes K a i t ; by calibration of the total electron 
transport rate estimated from fluorescence (/ raw ) with measure- 
ments of A+R d under non-photorespiratory conditions — where 
V. d \ t and V 0 are assumed to be absent — a calibrated electron 
transport rate (/ cal ) can be obtained. Under photorespiratory 
conditions J ca \ then represents the sum of V c and V 0 , such that: 

V 0 = /cal/4 -V c = ^(/cal/4 - {A + R d )) (5) 

From a theoretical perspective C c is relatively well defined, but 
see Parkhurst (1994) and Evans (2009) for issues with describing 
C0 2 fluxes or fluorescence with an average number representing 
different depths in the leaf. However, it is the practical estimation 
of / ca i and F a it that remains controversial and which represents 
a potential source of error in the calculation of g m . To obtain an 
accurate value for / cal , the raw measurements of chlorophyll 
fluorescence (/ raw ) must be calibrated and in doing so account for 
V a \ t under the experimental conditions as follows. Fluorescence of 
PSII provides an initial estimate of total electron flux through the 
electron transport chain: 

/ raw = O.425PPFD0psii (6) 

where 0.425 is the product of 0.85, the standard assumed value for 
leaf absorptance (a), and 0.5, the standard fraction of quanta 
absorbed by PSII relative to PSI (P), and <D PS n the quantum 
efficiency of PSII measured from fluorescence [®psn=(F m ' -F S ')I 
F m ']. If measured values for leaf absorptance are available, the 
assumed oc, and the estimate for / raw , can be improved. However 
the calibration procedures described below are often used to 
estimate a value for otp and therefore ot is not typically necessary. 
/ raw then can be related to /A+R d under appropriate non- 
photorespiratory conditions — normally at 1% 0 2 — where V 0 is 
negligible. From this relationship, the empirical values for a(3 can 
be found and thus provide a calibrated estimate of total electron 
flux (/ ca i)- Under non-photorespiratory conditions Equations 
5 and 6 become: 

/cai = 4(A + Rd) = mJ raw + c = m0.425PPFD(Dpsii + c (7) 

assuming a linear relationship. Thus the corrected value for otP is 
mxO.425. Alternatively, this equation is often converted from 



electron transport rates to quantum efficiencies by solving for 
^psii? preferably when no intercept is present: 

A + R d m0.425^ 
<J)C °2 = nnr7 ,^ = — ~ A — ®psii = m Opsn (8) 

where OC0 2 is the quantum efficiency of photosynthesis ([A+R d ]l 
PPFD), m' is the slope of the efficiency relationship, and the 
calibrated value for ocP is Am'. In practice, either of these 
relationships (Equation 7 or 8) are used for the calibration of / raw , 
with the fitted slopes providing an estimate of the value of ocP for 
the calibration conditions. The intercept is usually assumed to be 
zero. Alternatively, the presence of a non-zero j-intercept can be 
tested: if present, it represents alternative electron transport at the 
photo-compensation point. 

This calibration procedure is based upon the assumptions that: 

(i) ot and P are constant across the range of experimental variation; 

(ii) it is possible to estimate alternative electron transport as 
a constant proportion of total electron flux estimated as the 
intercept of the relationship; and (iii) the non-photorespiratory 
measurement conditions do not alter alternative electron transport 
relative to the experimental conditions. If either ocP or alternative 
electron transport vary with the environmental condition used to 
create the relationship (light or C0 2 ), non-linearities should be 
present in the relationship. An alternative is then to fit a non-linear 
function to the calibration data, such as the following linear- 
sigmoidal function: 

/cal = /A+R d = /raw ~ C ~ a /{I + exp[-(/ raw - b) / d]} (9) 



Analysis of sensitivity of 'variable J' g m magnitude to calibration 
scenarios 

To test whether calibration variants have significant effects on the 
calculation of g m , a sensitivity analysis was performed. Apparent 
shifts in g m due to changing the calibrations were calculated as 
follows for gas exchange measurements made on six leaves under 
ambient conditions (400 umol mol -1 C0 2 and a PPFD of 
400 umol m -2 s _1 ). (1) Standard calibration using assumed oc and 
P values (0.85 and 0.5) as is often used for the variable / method, 
with the following variants: (la) the mean R d and T* values 
measured using the Laisk method with six replicates; (lb) R d plus 
and (lc) R d minus the 95% CI of the mean; (Id) T* plus and (le) 
r*minus the 95% confidence interval of the mean. (2) Standard 
calibration using a measured oc (0.831) and assumed P value (0.5). 
(3) Calibrations fit to light response data measured under non- 
photorespiratory conditions at ambient C0 2 : (3a) using a linear fit, 
passing through the origin on the efficiency plot, but only using data 
points below a OC0 2 of 0.05 as suggested by Seaton and Walker 
(1990) and (3b) a linear-sigmoidal fit to the combined light response 
data on the rate plot. (4) Calibrations fit to the C0 2 response data 
measured at 400 umol m -2 s _1 PPFD and under non-photorespir- 
atory conditions, using the linear-sigmoidal function. Fitted param- 
eters for the calibration functions are provided in the Results. The 
non-linearity of the calibrations was assessed by comparing the 
Akaike Information Criterion (AIC) values between linear-sigmoi- 
dal fits and linear fits, where fits with the lowest AIC values 
have greatest support with model complexity taken into account 
(Burnham and Anderson, 2004). The R statistical program was used 
for these analyses (R_Development_Core_Team, 2010). An appar- 
ent value for g m was calculated for each of the six replicate leaves 
for all of the scenarios or parameter changes described above. 

Cross-validation of 'variable J' g m with g m estimated from curve- 
fitting procedures 

Values of mesophyll conductance to C0 2 were measured for an 
additional 10 C0 2 response curves using the same apparatus, 
corrections, and measurement conditions as detailed above. Added 
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to the five initial C0 2 response curves measured at 21% 0 2 , these 
provided a total of 15 curves, with an average of 14 C0 2 levels per 
curve. The measurements for the C0 2 response curves were made 
simultaneously with the fluorescence measurements, by using the 
2 cm 2 LI-COR fluorescence chamber, a necessary compromise, as 
the goal of this experiment was cross-validation between the 
variable / and curve-fitting methods. Ideally, measurements for 
use in curve fitting should be made using larger leaf areas (Warren, 
2006; Pons et al, 2009). 

The Exhaustive Dual Optimization (EDO) curve-fitting tech- 
nique of Gu et al (2010), as implemented on the LeafWeb website, 
was employed for this analysis in cognizance of the curve-fitting 
parameterization issues raised in that paper. The technique is 
based upon the principle that the photosynthetic C0 2 response 
curve can be represented by the minimum of a combination of 
three equations (Equation 10). These equations are non-rectangular 
hyperbolas that explicitly account for a non-infinite g m , and are 
based upon the original Farquhar-von Caemmerer-Berry-type 
photosynthetic functions: 

A =mm(A c ,Aj,A T pu) (10) 

The EDO approach uses functions for these three processes 
approximately similar to past curve-fitting approaches (Ethier 
et al, 2006; Warren, 2006), but applies these by assessing the 
possibility that any C0 2 response curve point could be limited by 
any of the three processes, with some constraints. The technique 
then exhaustively searches for parameter estimates for all of these 
possible limitation states, and selects the optimal fit as the fit with 
a minimum of a cost function consistent with the form of Equation 
10. For the EDO analysis, five parameters (g m , V c , max , J, R d , and 
the rate of triose phosphate utilization) were fit. To test for 
reliability of the parameter estimate for g m , the first and second 
derivatives of the cost function with respect to g m were tested to be 
zero and non-zero, respectively (Gu et al, 2010). 

To enable comparison between the two methods, g m values for 
the curve-fitting procedure (representing the entire C0 2 response 
curve) were compared with variable J g m values measured at 
ambient C0 2 (ambient C0 2 point on the C0 2 response curve), or 
the interpolated g m value for a Q of 600 umol mol -1 (interpolated 
as the point at a Q of 600 umol mol -1 on the line connecting the 
measured g m and Q value greater and less than 600 umol mol -1 ). 

Sensitivity analysis of 'variable J' g m to C0 2 variation 

The sensitivity of g m response to variation in C0 2 to errors in the 
estimation of R d or J cal was assessed by introducing a constant offset 
into a Farquhar-von Caemmerer-Berry-type photosynthetic model 
that held g m constant. A region of the modelled g m to Q response 
curve was defined from a Q of 200-500 umol mol -1 corresponding 
to values for which photorespiration should not be greatly inhibited, 
and thus measurements of V 0 would be relatively accurate. A linear 
slope was fit to these data, and non-zero slopes were recorded in 
response to introducing positive or negative biases in R d or / cal . 

The photosynthetic modelling was conducted using inputs of 
varying C c and constant values of g m (0.3 mol m - s -1 ), 
R d (0.42 umol m" 2 s" 1 ), r* (36.0 umol mol" 1 ), F cmax 
(70 umol m -2 s _1 ), and / (108 umol m -2 s -1 ). The values for these 
parameters were chosen to represent approximately a measured C0 2 
response curve for P. trichocarpa. The model calculations are 
provided online as a spreadsheet (Supplementary Spreadsheet SI 
available at JXB online). K c and K 0 values and the standard 
Farquhar-von Caemmerer-Berry equations were taken from von 
Caemmerer (2000), and reference to them is provided in the 
spreadsheet. From these inputs, V c ^ m [ n was calculated from the 
limiting process, namely the minimum of V c>c and V c j (the Rubisco 
and RuBP regeneration-limited carboxylation rates), and the photo- 
respiration rate, V Q , calculated from V c ^ n . The total RuBP 
regeneration rate, Jtotai/4, was calculated as the sum of K c?min and V 0 , 
assuming no alternative electron transport sinks and strict linear 



correspondence to fluorescence. Thus / to tai/4 provides a value for J C J 
4, as per the variable / method. Q was calculated using Fick's law 
and the calculated value for A. This model explicitly held g m 
constant; a bias was then introduced into the assumed value of R d 
or modelled value for J C JA (Jtotai/4), and from these new values 
a 'biased' estimate of g m was obtained using the formulae associated 
with the variable J method; the reverse of the initial calculations. 
Note that the formulae used in the variable / method given above 
(Equations 1^4) are algebraically the same as those used in the 
photosynthetic model just outlined. Thus the only difference 
between the biased estimate of g m and the value for g m when held 
constant is the introduction of a constant error in R d or J ca J4. 

Results 

Calibration of fluorescence estimates of electron 
transport rate with gas exchange 

Data used to calibrate fluorescence with gas exchange can be 
expressed either as quantum efficiency plots or as rate plots 
using C0 2 or electron equivalent units. As each has its 
advantages, the same light or C0 2 response curve data were 
compared on both plots (Fig. 1). The calibration relation- 
ships relating photosynthetic rates in electron-equivalent 
units (Ja+rJ to uncalibrated fluorescence estimates of 
electron transport (/ raw ) under non-photorespiratory condi- 
tions were non-linear when measured across a broad range of 
light or C0 2 conditions (Fig. 1). Light response curves 
demonstrated three phases of non-linearity: a subtle increase 
in / raw relative to J A +Rd at low light, a large shift towards 
increased «/ raw , but not J a +R d , at intermediate light, and in 
some responses a return to the one-to-one line at the highest 
light levels (Fig. 1A). For the same data plotted as efficiency 
plots (note the reverse in direction representing increasing 
light), the same shifts resulted in curvature towards greater 
quantum efficiency of net photosynthesis (®C0 2 ) at low light 
(Fig. 1C). This method of plotting the same data emphasizes 
the second curvature towards greater PSII efficiency (O PS h) 
at very low light (—100 umol m~ 2 s -1 ). 

Due to the curvature of the light response data from low to 
high light, two calibration relationships were fit. In the first 
method, a linear calibration was fit to each replicate light 
response on the efficiency plot forcing each line to pass 
through the origin (intercepts were not significantly different 
from the origin over this range of PPFD) and using data 
below a OC0 2 of 0.05 (here an average PPFD of >500 umol 
m~ 2 s -1 ), consistent with the suggestions of Seaton and 
Walker (1990) and resulting in a(3=0.383. If the calibration 
was done using an assumed value for ot of 0.85, as is common, 
the calibrated P value would be 0.451 rather than 0.5. The 
measured value of a was 0.831, resulting in an estimate of 
0.461 for the P value for higher light intensities. For 
comparative purposes, the one-to-one line was considered as 
the standard 'calibration' (ocP=0.425), as it is common to 
assume this value for aP with no further calibration. In the 
second calibration, a linear-sigmoidal curve was fit to all of 
the light response replicates simultaneously for the rate plot 
(Equation 9, a linear-sigmoidal fit: <z=ll.l, b=99.9, c=1.87, 
d=S.3l, adjusted i? 2 =0.974, AIC value=597.47). A linear fit 



418 I Gilbert et al. 




0 50 100 150 200 0 50 100 150 200 

J raw (umol e" m" 2 s- 1 ) J raw (jumol e" m" 2 S" 1 ) 




Fig. 1. Calibration plots for the rate of photosynthesis versus electron transport estimated from fluorescence (A and B), or for the 
photosynthetic versus fluorescence quantum efficiencies (C and D), measured under non-photorespiratory conditions (1% 0 2 ) using light 
response curves (A and C) or C0 2 response curves with 400 junol m -2 s~ 1 or 1000 jamol m -2 s~ 1 PPFD (grey and white symbols) 
(B and D). J raw was calculated using standard parameters (J raw =oepxPPFDx<Dp SN , ot=0.85, p=0.5, thus otP=0.425). Three lines are 
shown: the line where J ca \=Jmw (solid line; all panels), the average linear fit for nine light responses on the efficiency plot for points below 
a <DC0 2 of 0.05 (dashed line; C) and the average fitted linear-sigmoidal curve fit to the data of the nine light or five C0 2 responses (A and 
B). Different symbols represent measurements on different leaves. In C, one representative response curve is highlighted in black to 
illustrate regions of concave curvature at low cDC0 2 s and a final increase in cd psn at high <bC0 2 (low light). Arrows demonstrate the 
direction in which light or C0 2 increases on the different calibration plots. 



to the same data resulted in a marginally higher AIC value 
(598.81) and similar adjusted R 2 (0.974). 

It is common to use light response curves to calibrate the 
variable J method and then make use of them in studies 
using other experimental stimuli, for example variation in 
CO2. The carbon dioxide response curves measured at low 
0 2 and 400 umol m~ 2 s -1 PPFD did not resemble the light 
response curves at lower C0 2 concentrations, or at higher 
light (Fig. IB). A large shift was observed in the opposite 
direction to the light response curves, consistent with an 
increase in alternative electron sinks as may be expected by 
the lack of photosynthetic ability to use reductant under 



conditions of high light, low C0 2 , and low 0 2 . As the high 
C0 2 points were measured after the low C0 2 points, and 
showed high efficiency nearing the one-to-one line, photo- 
inhibition was not apparent (also F v '/F m ' returned to pre- 
low C0 2 exposure levels). The relationship on the efficiency 
plot was not linear; therefore, the C0 2 response curves were 
only calibrated using a linear-sigmoidal function on the rate 
plots (Equation 9, a linear-sigmoidal fit: a=22.8, &=95.3, 
c= -0.51, d= -8.6, adjusted i? 2 =0.958, AIC value=357.7). 
A linear fit to the same data resulted in a considerably 
higher AIC value (367.6) and lower adjusted R 2 (0.949), the 
difference between AIC values of —10 signifying that the 
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linear-sigmoidal fit had more support than the linear fit, ^ 
despite taking into account the extra parameters in the '<o 
linear-sigmoidal model (Burnham and Anderson, 2004). At 
a higher PPFD of 1000 umol m~ 2 s -1 , C0 2 responses cvj 
showed greater deviation from the one-to-one line, with q 
—40 umol e~ m~ 2 s -1 of apparent electron transport for -g 
little assimilation at the lowest C0 2 levels. To the authors' E 
knowledge, robust data have not been presented to validate ' 
that the calibration relationship is the same between ^ 
conditions of varying light and C0 2 , apart from Hassiotou 1 
et al (2009) whose results largely confirm those in Fig. 1. 
Most studies use a single saturating flash to measure F m \ *v 
potentially introducing additional non-linear effects with 
changing light (Markgraf and Berry, 1990; Earl and Ennahli, 
2004). However, measurements of varying light demonstrated 
that using multiple saturating flashes rather than a single 
flash did not linearize the calibration functions, and rather 
the size of the discrepancy between /A+R d and / raw was 
slightly enhanced at high PPFDs (F m 9 increases while F s 
remains constant). As the key C0 2 calibrations were done at 
moderate PPFD (400 umol m~ 2 s -1 ), these were not affected. 



Size of calibration biases on J ca/ and g m 

From the data shown in Fig. 1, five types of calibrations 
were performed: the three linear or linear-sigmoidal func- 
tions, the standard assumed calibration parameters from 
the literature, and the measured leaf absorbtance and 
assumed (3. The magnitude of the errors in J cal on the rate 
plots, and particularly the efficiency plots, is both difficult 
to visualize and hard to relate to the magnitude of the 
measured quantities. Therefore, data were expressed as the 
residuals for the rate relationship, rescaled to units of C0 2 
uptake, and plotted against the PPFD or C0 2 used to 
generate the points (Fig. 2 A, B). The residuals were 
calculated as J ca \/4-A, where J ca \ is the electron transport 
rate calibrated using one of the five types of calibration. 

The linear-sigmoidal calibrations applied to the same 
environmental variation to which they were fit produced the 
smallest residuals, and did not have any systematic errors 
across a broad range of light or C0 2 , possibly apart from 
2000 jimol m" 2 s" 1 PPFD (Fig. 2A). The linear higher light 
calibration produced few residuals at high light, but 
consistently underestimated A by ~1 jimol C0 2 m~ 2 s -1 at 
low light. The standard calibration, assuming a(3=0.425, 
performed poorly, with significant overestimates of A of up 
to 4 jimol C0 2 m~ 2 s -1 under all but low light conditions. 
The standard calibration was even worse for low C0 2 
conditions (Fig. 2B), resulting in residuals as high as 6 jimol 
C0 2 m~ 2 s -1 . The calibration with assumed value for (3 and 
measured a had a similar pattern to the standard calibration 
although the residuals were improved. 

The variation between these calibration curves resulted in 
large differences in apparent g m values when applied to 
measured photosynthetic data for ambient C0 2 and moder- 
ate light (Table 1). Values ranged by 104% from a minimum 
of the standard assumed calibration to that of the linear 
calibration. Linear-sigmoidal fits to light or C0 2 response 
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Fig. 2. The average residuals of the calibrated rate of electron 
transport (J ca! ), rescaled to units of C0 2 uptake, relative to the 
observed photosynthetic rate for light (A) or C0 2 response curves 
(B) measured under non-photorespiratory conditions. To allow for 
possible trends in R d , residuals were calculated as J ca /4-A Points 
represent the mean and standard errors for five or more replicate 
light curves (the same data as in Fig. 1 A and C) or five C0 2 
response curves (the same data as in Fig. 1 B and D). 

curves were intermediate. This variation in apparent g m due 
to the underlying calibration was larger than variation in g m 
caused by changes in ^ or when adjusted by the 95% 
CIs of the mean values (Table 1). 

Correspondence between 'variable J' g m and curve- 
fitting g m 

The g m values calculated by the variable J and EDO 
approach curve-fitting method were most highly correlated, 
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Fig. 3. Cross-validation of g m values calculated from 15 C0 2 
response curves using three alternative calibrations for the variable 
J method and applied to the ambient C0 2 measurement on the 
curve, and g m calculated from the Exhaustive Dual Optimization 
(EDO) approach for fitting Farquhar-von Caemmerer-Berry models 
of Gu et al. (2010). Experimental conditions were: PPFD, 400 jimol 
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s" 1 ; T h 24.9±0.8 °C; VPD, 1.47±0.49 kPa. 



and the points nearest the one-to-one line for the C0 2 
linear-sigmoidal calibration (Fig. 3). The assumed and 
linear-light calibrations resulted in correlations between the 
variable J g m and the curve-fitting g m values, but resulted in 
greater deviation from the one-to-one relationship. In 
addition, the linear-light calibration resulted in a negative 
value of g m , which was removed from the analysis. Two g m 
values were removed from all analyses due to the EDO 
curve-fitting analysis providing high g m values (>0.5 mol 
m~ 2 s -1 ), and this was consistent with a zero second 
derivative of the EDO cost function (the condition under 
which the parameter estimate is not reliable). For this 
analysis the variable J g m estimate was limited to measure- 
ments made at ambient C0 2 , while the curve-fitting estimate 
used the entire C0 2 response curve data. When the variable 
J g m value representing a Q of 600 umol mol -1 was plotted 
against the curve-fit g m value, R 2 s were reduced and the 
variable J g m value was an underestimate for all calibra- 
tions. For the linear-light calibration a number of variable J 
g m estimates at high Q were negative. 

Response ofg m to C0 2 

The response of g m to C0 2 (detrended for stomatal 
conductance changes by using Q) was highly variable when 
the five calibration protocols were compared (Fig. 4). In one 
of the five replicates (Fig. 4A), all five calibrations produced 
values for g m in the range of past reports (Niinemets et al, 
2009); in the other four replicates the linear and linear- 
sigmoidal fit to the light response calibration resulted in 
negative or large (>1) values for g m at C0 2 levels higher than 
ambient (one representative replicate is shown in Fig. 4B). 
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Fig. 4. The observed response of g m to C-, illustrating either 
qualitative agreement amongst the four calibration options (repre- 
sentative of only one of five replicate C0 2 response curves; A) or 
marked disagreement between calibration methods (a representa- 
tive C0 2 response curve for four of five replicates; B). Note that in 
B two of the calibration options result in negative values of g m due 
to overestimates of C c (apparent C c >Q). These occur as g m =A/ 
(C-C c ), thus underestimates of J ca! at high C-, due to the different 
calibrations lead to C c approaching C„ the denominator of the 
equation is small leading to g m approaching infinity, and when C c 
becomes higher than C„ g m instantly becomes negative. Experi- 
mental conditions were the same as in Fig. 3. 

As only the standard calibration constants and the 
linear-sigmoidal fit to the C0 2 response calibration gave 
reasonable values for g m for all replicates, these two 
calibration protocols were investigated in greater detail. 
Using either calibration, g m showed strong shifts, increasing 
from the lowest Q values, remaining stable or slowly 
decreasing at ambient CO2 values, and decreasing strongly 
at high Qs (Fig. 5A, B). However, the Harley et al (1992) 
criterion was violated for almost all points at high Q. 
Nevertheless, the points that satisfy the Harley criterion 
(Harley et al, 1992) demonstrate a consistent negative 
response of g m to Q. 
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Fig. 5. Observed g m response to C-, for five replicate C0 2 
response curve using the standard calibration constants (A) and 
the linear-sigmoidal fit to the C0 2 response calibration (B). Points 
that satisfy the Harley criterion (filled symbols) and points that had 
a Harley criterion of <10 or >50 (open symbols) are distinguished. 
Replicate curves are shown with the same symbols. Experimental 
conditions were the same as in Fig. 3. 

The potential for calibration biases to affect the calcula- 
tion of g m can be illustrated by plotting Q versus the 
parameters used to estimate C c and g m (Fig. 6). As Q 
increases, V 0 tends towards zero due to competitive in- 
hibition of photorespiration by C0 2 . As a result, C c 
— proportional to the VJV 0 ratio — is increasingly vulnera- 
ble to biases at high Q. Specifically, as V Q decreases at 
higher C0 2 s, any errors in its estimation lead to an inflated 
C c (calculated from the V C IV Q ratio), as C c tends towards Q, 
g m values [calculated from A/(C-C C )] rapidly become large. 
Once C c , estimated with slight errors, is the same or larger 
than Q, g m becomes infinite or negative. This explains the 
variability, high and negative values of g m in Fig. 4. 
Furthermore, as V 0 tends towards zero, biases in J ca \/4 due 
to alternative electron transport sinks (or changes in R d , a, 
or (3) become increasingly important. In other words, small 
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Fig, 6, Rates of the standard fluorescence estimate of electron 
transport (J raw ), the calibrated rate of electron transport expressed 
in C0 2 equivalents, carboxylation (V c ), apparent photorespiration 
rate (V Q apparent =J ca /4-\/ c ), and photorespiration calculated 
assuming C C =C; (V Q at infinite g m ), for a measured response to 
C0 2 . Note the small shift in V Q necessary to result in an infinite g m 
at high C0 2 . Points that satisfy the Harley criterion (filled symbols) 
and those that did not (open symbols) are distinguished. Experi- 
mental conditions were the same as in Fig. 3. 

errors in the estimation of J ca \ have increasing impact on the 
estimation of g m at high Q, as V Q becomes small and the 
error to the V Q ratio increases. 

Sensitivity of response of g m to C0 2 

The sensitivity of g m to calibration biases was further 
investigated by introducing small biases in J ca \/4 into 
a photosynthetic model. The lines in Fig. 7A represent the 
apparent g m response for simulated data for which g m was 
held constant, but for which small systematic errors were 
introduced into the value for J ca \/4 and g m then calculated 
from the biased data. Any overestimation of J ca i/4 results in 
a lower g m and an apparent negative relationship with 
increasing Q (Fig. 7A). In contrast, underestimating J ca \/4 
results in a larger apparent g m . The presence of a positive or 
negative relationship between g m and Q was a function of 
the small constant biases added to J C J4 (Fig. 7B). If, in the 
photosynthetic model, g m is assumed to be constant with 
C0 2 , then the residuals in Fig. 7B demonstrate that 
previously used calibration relationships would consistently 
result in apparent negative g m responses to CO2, while the 
linear-sigmoidal C0 2 response calibration would result in 
both negative and positive relationships. If this assumption 
is true, the sensitivity analysis demonstrates that the bias in 
J ca i/4 necessary to result in an artefactual effect of Q on g m 
is small (<0.5 umol C0 2 m~ 2 s -1 ) relative to the residuals 
typically observed in the calibration relationship (—0.5- 
6 umol C0 2 m" 2 s" 1 ; Fig. 2A or B and Fig. 7B). If, 
however, g m is truly not constant, the observed slope of the 
response of g m to C0 2 would still be sensitive to errors, and 
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Fig. 7. The sensitivity of the g m to Q response to constant errors 
added to J ca /4 for a photosynthesis model that had a constant g m 
(A) and the sensitivity of the slope of the g m to C-, response for the 
same errors in the photosynthesis model (B). Average and 
standard deviations of the residuals [J c J4-(A+R d )} of three 
calibration options are shown in B. The slope of the g m to Q 
response was defined as the linear fit to the modelled data within 
a range of Qs that would satisfy the criterion of Harley et al. (1 992) 
(200-500 umol mol -1 ). The residuals were calculated across the 
same range of C,. Numbers in A represent the overestimate (error) 
added to J ca /4 where zero error (and thus the true modelled 
relationship) had no relationship between g m and C h Vertical lines 
in A represent the C-, at which g m estimates shifted from high to 
negative values (see text for explanation). 

if errors were large enough would result in transitions from 
negative to positive or vice versa. The model is included as 
a spreadsheet in the Supplementary data at JXB online. 



Discussion 

Can the variable J method measure the response of g m 
to C0 2 ? 

The nature of the observed response of g m to C0 2 is highly 
sensitive to biases in the estimation of the calibrated total 



photosynthetic electron transport rate (J ca \). This was 
demonstrated using a sensitivity analysis of standard 
photosynthetic equations, to which a systematic bias was 
added. For example, the sensitivity is such that there is an 



apparent 23% decrease in g m over a 300 umol mol~ range 
of Q when an —2% (—0.5 umol C0 2 m~ 2 s -1 ) overestimate 
of Jcai/4 is included in the photosynthetic model, despite the 
modelled g m remaining constant (Fig. 7A). As the true 
modelled relationship was on the steepest portion of the 
sensitivity analysis (Fig. 7B), this demonstrates that if g m is 
indeed constant, then any bias in J ca \/4 will lead to 
artefactual positive or negative relationships of g m to Q. If 
g m is dynamic, varying with C0 2 , the point of greatest 
sensitivity will drift, but the overall pattern of sensitivity 
demonstrated will remain. In this case, an observed re- 
lationship may represent a true response, but the slope will 
be sensitive to measurement errors and calibration biases. 
This sensitivity analysis provides similar results to those 
which Harley et al (1992) presented in their fig. 6, and those 
which Hassiotou et al (2009) presented in their supplemen- 
tary material. Indeed, Harley et al (1992) note that: 'In all 
cases, the sensitivity to errors was relatively low between 
100 and 300 ubar Q, but outside this range the sensitivity 
was so great that the results could become unreliable.' 
Despite these earlier cautions, subsequent researchers have 
continued to use this approach over a broad range of 
conditions. It is important to note that these considerations 
are applicable to any environmental variation that may 
affect photorespiration: C0 2 , temperature, light, stomatal 
closure, etc. For instance, a similar analysis could be done 
for the relationship of g m to PPFD, in which case the 
relationship would be sensitive to errors at PPFDs below 
light saturation where the errors become significant relative 
to photorespiration. Thus it is also the relationship of g m to 
light that is sensitive to errors when using the variable J 
method, although at saturating light intensities the presence 
of high rates of photorespiration leads to less sensitive 
estimates of the relationship of g m to light. 

The residual variation in the different calibration rela- 
tionships — which is a determinant of the error in J C J 
4 — was up to 5 umol C0 2 m~ 2 s -1 at the extreme of using 
standard calibration constants, and about ±1 umol C0 2 
m~ 2 s -1 when calibrated using a linear-sigmoidal function 
on C0 2 response data (Fig. 2A, B). Thus, the magnitude of 
the errors in the calibrations was similar to, or considerably 
larger than, the error necessary to affect whether there is an 
apparent response of g m to Q (Fig. 7B). It is broadly true 
then, given the large errors in our best estimates of J CSL \/4, 
and the sensitivity of the g m to Q relationship to any error, 
that it is difficult to measure the response of g m to Q using 
the variable J method. That is, with the high overestimates 
of ,/ ca i demonstrated for standard calibration methods over 
a moderate C0 2 range (Fig. 7B), the variable J method is 
likely to produce steeper relationships between g m and C0 2 
than actually exist. 

Variable J g m and partially independent g m values from 
the EDO curve-fitting approach corresponded well when 
the variable / technique was limited to use under ambient 
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C0 2 and with the non-linear calibrations reported here (Fig. 
3). These results are consistent with the sensitivity analysis 
performed earlier (Table 1). That is, relative to the EDO 
curve-fitting g m values the linear-light calibration causes 
overestimates in variable J g m , the C0 2 linear-sigmoidal 
calibration results in approximate correspondence, and the 
assumed calibration results in underestimates. This suggests 
that when appropriately calibrated the variable J method 
has value for studies comparing species, using unstressed 
plants at moderate light and ambient C0 2 , but should not 
be applied across a range of environmental conditions 
under which photorespiration is likely to vary. However, 
these non-linear calibrations remain empirical and do not 
address the implication that a non-linear response repre- 
sents an unaccounted for fundamental change in photosyn- 
thetic functioning. 

Why are the calibrations non-linear? 

It is difficult to provide a retrospective review of whether 
the non-linearities in the calibration relationships observed 
here are present in the g m literature. For example, a litera- 
ture review of 56 experimental studies of g m , published since 
1992, found that 66% of these use the variable J method, 
and 44% use it as a sole technique. Of these, few studies 
provided calibration data, and if this was done even fewer 
calibrated the variable J method using environmental 
variation appropriate for the experiment at hand. Fewer 
performed brief sensitivity analyses, and finally no study 
attempted to calibrate the technique using non-linear 
functions. However, many of the non-linear effects de- 
scribed here have been previously described by Seaton and 
Walker (1990) and Oquist and Chow (1992). There are also 
indications of non-linearities in the calibration relationships 
used to calculate g m or C c (Warren, 2006; Galle et al, 2009; 
Hassiotou et al, 2009; Loreto et al, 2009). 

Seaton and Walker (1990) and Oquist and Chow (1992) 
demonstrated large non-linearities in light response curves, 
plotted as efficiency plots, measured under non-photores- 
piratory, saturating C0 2 conditions and with oxygen 
electrodes. These curved relationships on efficiency plots 
are consistent with the sigmoidal patterns found on the rate 
plots, but there are clear differences in weighting of points 
between the plots. The reasons for the non-linearities are 
discussed by Oquist and Chow (1992) and include: (i) 
changing connectivity of PSII units, leading to more cycling 
of electrons between chlorophylls; (ii) at low light, mito- 
chondrial respiration (R d ) may increase, but in the variable 
/ calculations Rd is assumed to be constant and a single 
value usually estimated for all conditions from Laisk 
response curves; (iii) fluorescence parameters may be 
estimated from slightly shallower populations of chloro- 
plasts than those that fix C0 2 , and the contributions of 
these populations of chloroplasts would change with light 
intensity (Warren, 2006; Evans, 2009). On the rate plots, 
possible alternative electron sinks are highlighted, resulting 
in non-linear shifts in the calibration relationship, and may 
represent little (Ruuska et al, 2000), or up to 24% of the 



total electron flux (Haupt-Herting and Fock, 2002). Two 
main processes are thought to account for alternative 
electron sinks (von Caemmerer, 2000), each accounting for 
up to 10% of total electron flux: reductant provided to 
nitrate assimilation (Rachmilevitch et al, 2004) and the 
Mehler reaction (Haupt-Herting and Fock, 2002). 

These effects are highlighted when comparing light and 
C0 2 responses measured under non-photorespiratory con- 
ditions (Fig. 1A, B). A priori this must be expected, as at 
low C0 2 , and particularly at high light, there is a limitation 
on reductant use, but high reductant supply that will result 
in large alterations of PSII heat dissipation and may result 
in up-regulation of alternative dissipative energy sinks, such 
as the Mehler reaction (Neubauer and Yamamoto, 1992). 
The quantitative effects of alternative energy sinks remain 
unclear (Ruuska et al, 2000); however, it is noted that 
relative to the errors (—2% of /totai/4) necessary to cause 
apparent changes in g m , estimates of alternative electron 
sinks are large and therefore vital to account for. 

Finally, it is not clear whether alternative electron sinks 
are changed when shifting from ambient to low 0 2 as 
required for the calibration curves (Pons et al, 2009). For 
instance, at high C0 2 the calibration curve was closer to the 
one-to-one line than for high light points (Fig. 1). This may 
imply that alternative electron transport sinks are affected 
by the capacity of photosynthesis to dissipate absorbed light 
energy, or are directly affected by C0 2 . Considerable shifts 
in nitrate assimilation with age, C0 2 , and oxygen concen- 
tration occur in Arabidopsis, using an equivalent electron 
flux up to 10% of the photosynthetic rate (Rachmilevitch 
et al, 2004), and thus represent evidence of alternative 
electron transport shifts that could occur during the 
calibration procedure. If this is generally the case, it would 
be challenging to find conditions under which the variable 
/ method can be calibrated. Indeed, the fitted a or (3 
parameters for a non-linear calibration function cannot 
then be interpreted as physical constants as the non- 
linearity implies that they change with environmental 
conditions, or that they include alternative electron trans- 
port sinks. It appears that much work remains to be done, 
using independent methods, to understand the implications 
of the photosynthetic changes that occur when producing 
calibration relationships for the estimation of g m and C c 
using the variable J method. 

How should the variable J method be used? 

The variable J method appears difficult to validate under 
circumstances of varying photorespiration due to the 
extreme sensitivity of g m under conditions of low photores- 
piration. However, the method when calibrated taking non- 
linearities into account did improve the estimates of g m 
under ambient C0 2 relative to the EDO curve-fitting 
approach. Thus if the variable J method is to be used for 
comparing species (and not environmental variation) the 
following are imperative: (i) a calibration is done with 
conditions that match the experimental conditions (not 
a light calibration versus C0 2 experiment); (ii) the 
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calibration (and experiment) is limited to the linear region, 
for example OCO 2 <0.05; Seaton and Walker (1990), or non- 
linear functions are used, and if not linearity should be 
explicitly tested; (iii) the calibration is fit using rate plots, 
not the efficiency plots that emphasize low photosynthetic 
rate points disproportionately; and (iv) a sensitivity analysis 
is done that asks what size biases in the estimation of J ca \, 
or variation in the values for R d and T*, are necessary to 
remove the observed effect or relationship, and are such 
errors plausible for the calibrations. Regardless of these 
improvements, the lack of knowledge of why the calibration 
response is curved, and whether alternative electron sinks 
are affected by changing 0 2 may preclude the use of the 
variable J method in most experiments. 



Conclusion 

The variable J method is sensitive to errors and must be 
used with caution in experiments where photorespiration 
varies. Nevertheless, none of the calibration or sensitivity 
scenarios tested here precludes an effect of any variable on 
g m \ thus g m may be dynamic rather than constitutive, but 
these results suggest that we cannot know the magnitude or 
nature of changes with certainty using this technique. It is 
suggested to limit use of the variable J method to 
comparing species under conditions of moderate light and 
ambient C0 2 with appropriate calibration, and not in 
experiments measuring responses to environmental factors 
that affect photorespiration. There is much research needed 
using independent methods to provide information on 
whether and how g m and alternative electron sinks respond 
to C0 2 , light, or 0 2 . The region at which g m measured using 
the variable J method starts declining with C0 2 (Flexas 
et al, 2007) or reduced light (personal observation) 
corresponds to the point where RuBP regeneration becomes 
limiting to photosynthesis. Although this may occur 
through a common mechanism related to RuBP regenera- 
tion, this effect is less apparent in g m measurements using 
carbon isotope discrimination (Flexas et al, 2007; Tazoe 
et al, 2009, 2011; Vrabl et al, 2009). The point where RuBP 
regeneration becomes limiting for both the light and C0 2 
response curves also corresponds to a decrease in photores- 
piration. Thus at this point the ratio of biases to photores- 
piration dramatically increases, causing artefacts to be 
introduced into the response of variable J g m to C0 2 or 
light, if subtle biases are present in the calibration or 
measurements. Thus, it is suggested that positive biases in 
the calibration procedure result in the variable J method 
overestimating the slope of the relationship between g m and 
C l — an explanation for the differences between studies using 
the variable J method and those using carbon dioxide 
discrimination. 

Supplementary data 

Supplementary data are available at JXB online. 



Spreadsheet SI. Sensitivity analysis of modelled photo- 
synthetic response to C0 2 , with g m held constant. 
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